Hubbard ladders in a magnetic field 
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The behavior of a two leg Hubbard ladder in the presence of a magnetic field is studied by means 
of Abelian bosonization. We predict the appearance of a new (doping dependent) plateau in the 
magnetization curve of a doped 2— leg spin ladder in a wide range of couplings. We also discuss the 
extension to N— leg Hubbard ladders. 
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In the past few years there has been intense research 
activity, both experimental and theoretical, on ladder 
systems These systems interpolate between Id and 
2d and, due to the strong quantum fluctuations enhanced 
by the low dimensionality, exhibit surprising and exotic 
behaviors, such as, for instance, the so-called even-odd 
effect. 

A feature of spin ladders which attracted considerable 
interest is the occurrence of plateaux in their magneti- 
zation curves. One of the main observations has been 
the rationality of the value of the magnetization at which 
plateaux can appear However, it has been recently 

pointed out that plateaux can also appear at irrational 
values of the magnetization. This happens for example 
in periodically modulated doped Hubbard chains J5[ , but 
this scenario is believed to be generic in doped systems. 
A novel and interesting feature of the plateaux predicted 
in B is the fact that their position depends on the filling 
and therefore can be moved down to low magnetization 
values by means of doping. This makes doping depen- 
dent plateaux potentially observable in experiments at 
low magnetic fields. Previous studies revealed doping 
dependent plateaux in other systems, such as the one- 
dimensional Kondo lattice model || and an integrable 
spin-S" doped t — J chain Q| . 

The purpose of this letter is to investigate the condi- 
tions of occurrence of magnetization plateaux in doped 
ladder systems. We have focussed on the case of cou- 
pled Hubbard chains, the spin ladders being particular 
limits of Hubbard ladders at half filling. Various mod- 
els of coupled chains have been thoroughly investigated 
in the last decade, but comparatively little attention has 
been devoted to the effects of a magnetic field on their 
properties. In [|| (see also ||) it has been proven that in 
a N-\eg Hubbard ladder gapless excitations exist at mo- 
mentum k = 2itNri'\.i if the commensurability conditions 
Nn^i e Z are not satisfied JhJ. This holds even if the 
low energy excitations cannot be identified as the usual 
separated charge and spin excitations. According to this 
result, magnetization plateaux at to ^ can only occur 



if at least one of the conditions Nn^i = |(n±m) e Z is 
satisfied, which considerably restricts the window to find 
new plateaux. In order to prove the presence of a plateau 
one has to show that there is a finite gap in the total spin 
sector, since the above requirement just provides a neces- 
sary condition for the occurrence of a non trivial to =/= 
plateau. We have found that when the conditions 



N 

Nn u = j(n±m) e Z 



(1) 



are simultaneously satisfied, or when only one of them is 
satisfied and the doping kept fixed, a plateau can indeed 
occur for some range of parameters. In the last situation, 
one has a doping dependent plateau which, for the spe- 
cific case of N — 2, is predicted to be present for a wide 
range of couplings. 
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FIG. 1. Schematic magnetization curve of a 2-leg Hubbard 
ladder at n = 1 (half-filling) and n 7^ 1. The width of the 
plateaux depends on t±/t, U/t,n. 

Our main result is schematically summarized in Fig. 1, 
which shows the expected magnetization curves of a 2-leg 
Hubbard ladder for both the doped and undoped cases. 
According to our analysis, apart from the known plateau 
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at m = 0, plateaux should also occur at m — ±(1 — n). 
These could be observed experimentally by doping any 
system described by a 2-leg ladder. Possible candidates 
are notably the spin-i ladder compounds SrCu2C>3 fTl| l 
and Cu2(C5Hi2N2)2Ci4 (which has been studied under 
high magnetic fields Jl2]| ) , or some quasi-one-dimensional 
spin-i organic system. These materials are often difficult 
to dope. Nevertheless, doping has been achieved for the 
ladder material Sro.4Ca13.6Cu24O41.s4 under high pres- 
sure |l3| (in the latter, superconductivity was also ob- 
served). Another class of interesting candidates are sin- 
gle wall carbon nanotubes Jl4| , whose low energy physics 
is also described by a model equivalent to two coupled 
spin- 1/2 fermionic chains away half filling pq |. 

Let us turn now to our bosonization analysis (see e.g. 
frill). We start from the 2-leg Hubbard ladder, with the 
following lattice Hamiltonian: 

H = H I + H II + \ C J2 n WJ + ~ X sJ2 S * S " 

X X 

X 

where Hi (i = /, II) is the usual Hubbard Hamiltonian: 

Hi = -t 5^( c 2+l,a4,a + H.C.) + <A<1 + 



T "r "x,l 



)-oE«T"<i)' ( 3 ) 



, c x are electron creation and annihilation operators 
on site x of chain i, a is the spin index, n x a = c l J a c xa . 
The index i = I, II can be considered either as a chain 
index or as a band index. Accordingly, the Hamiltonian 
(||) corresponds to different regimes of coupled chains of 
interacting electrons. These systems have been largely 
studied analytically, for different range of parameters 
0- |5[, and numerically ||-|§]. 

Eq. (0) describes two identical chains with standard 
Hubbard Hamiltonians Hi^it, Uq, /^o, h) with small in- 
trachain repulsion Uq, coupled by a direct hopping term: 



Hi, 



:t c 2 

X.OC X.Oi 



H.c.) , 



(4) 



where t± is the interchain hopping amplitude. To see 
this, one simply has to change variables to the bonding 
and anti-bonding basis: 



V2 



V2 



(5) 



Rewriting H\-\-Hi-\-Hi nt in this basis and identifying the 
band indices s and a with the indices / and II, one recov- 
ers the Hamiltonian (||) with fJ.1,11 = fi^f t±, U = Uq/2, 
4A C = — X s = — 2A = Uq. Thus we have to deal with two 
coupled inequivalent effective Hubbard chains. 



On the other hand, for Uq 3> t± and close to half filling, 
(namely, for strong repulsive interactions in the individ- 
ual chains), the effect of band splitting is suppressed due 
to the requirement of avoiding double on-site occupancy. 
In this situation it is more appropriate to consider cou- 
plings directly in terms of the original degrees of freedom 
of the two individual chains; then, / and iT have to be 
identified with the chain index, leading us to two equiv- 
alent chains coupled by density-density and spin-spin in- 
teractions, as in ^ but with A = and A CiS ~ t\/U. 

We are now ready to treat the two different situa- 
tions at once, by studying the Hamiltonian (|^), which 
contains the two limits described above. The m = 
case has been largely studied in the literature JTtJ- [p8[ . 
The main conclusions in that case are that the spin gap 
present at half filling survives upon doping, although 
smaller in magnitude, and the appearance of quasi-long- 
range superconducting pairing correlations. The main 
effect of the magnetic field on the single chain is to shift 
by opposite amounts the up and down Fermi momenta, 
whose difference is proportional to the magnetization, 



-fc. 



irm. As a first approximation, we can ne- 



glect other effects, such as the mixing of spin and charge 
degrees of freedom, and the dependence of the scaling 
dimensions of the operators on the magnetic field p9[ . 
Note that the calculations can also be handled straight- 
forwardly without these approximations, but equations 
become cumbersome (see for example p|). 

In order to write the bosonized expressions for the two 
cases for m / we have to define the Fermi momenta 
/c|'| 7 separately for each band (or chain) and linearize 
the dispersion relations around them. We assume that 
interactions do not change significantly the free Fermi 
momenta. In the process of bosonization we obtain oper- 
ators containing oscillating factors made of various com- 
binations of the k l a . In the case of inequivalent chains, 
most of them will be commensurate only if specific con- 
ditions depending on t± are satisfied. We will neglect in 
the following these non-generic situations and will take 
into account only operators that can be commensurate 
independently of the value of t±. Moreover, the operator 
responsible for the m = plateau is only commensurate 
for that value of m and hence it will not enter in the 
following discussion. 

i) Equivalent chains 

In this case = k^ 1 , I, II represent chain indices 



and we use the notation I = 1,11 
Hamiltonians Hi are written as 



2. The continuum 



(6) 



where V c , s ,Kc,s are the velocities and the effective Lut- 
tinger parameters for charge and spin degrees of freedom. 
Away from half-filling (n ^ 1) and for non zero magneti- 
zation J30|), the effective interaction Hamiltonian reads 



2 



Him = A c d x & a d x & c + -j d x *%& B + 

- Ai (— cos[\/47r$ ( :] cos[\/47r$7] - 
+ cos[2(fc{ + kf)x - VSn^} + 
+ cos[2(fc{ + k^)x - %/8^$+] \ + 
+ A 2 cos[%/47r$~] cos[\/47r67] 7 



(7) 



where $ l c/s = ±*j), and *± = ^ ($ x ± $ 2 ) ||. 

We have kept only the more relevant operators for the 
case of large U. 

The bosonic bilinear terms can be absorbed in the ki- 
netic part of the Hamiltonian by moving to the ± basis 
defined above. As a consequence, the velocities and more 
importantly the effective Luttinger parameters are renor- 
malized in the following way: 



K c -> Kf = K c {\ ± KKc/vc)- 1 ' 2 , 
K s - Kf = K s (l ± X.K.lAv,)- 1 /* 



(8) 
(9) 



From the above equations we see that in the charge sec- 
tor, the Luttinger parameter for the total charge field 
(symmetric combination) is slightly reduced with respect 
to K c whereas the one for the relative charge (anti- 
symmetric combination) is slightly increased (we assume 
A c > 0). Let us recall that K c decreases from 1 for U = 
to 1/2 for infinite repulsion, thus if+ < 1. We can repeat 
the analysis above for the spin sector and the scaling di- 
mensions of the perturbing operators in eq. ^ at zero 
loop are then given by: 



(10) 

(11) 
(12) 



cos[v47r$ c ] cos[V47T $ s ] -> K c + K 8 

K+ 



/8^$+ T ] 



id 



cos[v47r$~] cos[V47r87] -> K~ + 1/K~ 



Due to the double cosine terms in eq. (Q) we conclude that 
both spin and charge antisymmetric sectors are massive. 
Let us now consider the symmetric sectors, which are af- 
fected by the terms in the third and fourth lines of eq. 
( 7 ) . These operators are commensurate only if conditions 
(1) are satisfied. These appear to be exactly the same 
conditions we found in the dimerized Hubbard chain for 
the appearance of plateaux || . When the two conditions 
are simultaneously satisfied we expect both and 
to be massive. Indeed, both and are decreased 
with respect to their values in the absence of interchain 
coupling, and we can conclude that the operators, when 
commensurate, are relevant |32 ]. When instead only one 
of these conditions is satisfied, let us say the one with the 
+ sign, only the field gets a relevant cosine interac- 
tion. Then following the analysis in J5| we can conclude 
that a magnetization plateau occurs at m — 1 — n if the 
total density is kept fixed. 

ii) Inequivalent chains 



In this case, /, II are band indices. When t± > 2t, the 
bands are splitted and the situation is completely analo- 
gous to one found in the doped p-merized Hubbard chain 
we treated in ||. Following the same lines, we can prove 
that a plateau occurs when 2n^i € Z. Notice that in 
this particular regime the plateau at m = is no longer 
present. When t± ~ t, non trivial processes between 
bands are allowed and in this case the effective interac- 
tion Hamiltonian away from half-filling (and at m =/= 0) 
reads 

H int = Xcd^id^i 1 + ^ d^id^ 1 / + 

- Ai ( cos[2(fcf + fef )x - \/8^] + 
+ cos[2(fc{ + kY)x - V8^<I>+] ) + 
+ A 3 cos[\/47r07] cos[\Z47r0,7] . 



(13) 



In this expression, we have neglected all terms containing 
combinations of k l a de pen ding explicitly on t±. The last 
marginal term in eq. (03) is generated radiatively. The 
derivative terms lead to a renormalization of the Lut- 
tinger parameters as in eqs. (|[^). 

As in the strong U limit, when 2n^i £ Z, the sec- 
ond and third line of (|l3| ) become commensurate. We 
can now repeat the analysis concerning the relevance of 
these operators. Their dimensions are determined by the 
forward scattering terms which can be computed at first 
order in U -C t, t±. It is straightforward to see that K+ s 
is renormalizcd from 1 to < 1, which makes the 



term cos[v87r$| ] (respectively cos[V87r$| ]) marginally 
relevant. From this and the previous analysis, we can 
therefore conclude that both in weak and in strong cou- 
pling regime of the 2-leg Hubbard ladder, a magnetiza- 
tion plateau is expected when 2n-|- 1 | S Z, and the magne- 
tization curve at density n < 1 should have the schematic 
form depicted in Fig. 1. 

It is interesting to analyze the behaviour of the diago- 
nal components of the charge density wave (CDW) and 
superconducting pairing (SC) operators at a doping de- 
pendent plateau, as it occurs e.g. for m — 1 — n. In the 
case of equivalent chains the behaviour of the correlators 
of such order parameters depends on which operator, be- 
tween the first and the last double cosine terms in eq. (Q) 
dominates. When the first term dominates, one could 
find among diagonal SC and CDW operators p3j 



CDW,i 



i\/27re 



+ 



(14) 



algebraic decaying correlators. On the contrary, when 
it is the last one which dominates, then all correlators 
decay exponentially. One can argue that varying the val- 
ues of A s and the magnetization, both situations could 
be achieved, but this requires a more detailed analysis of 
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the Bethe Ansatz equations in the presence of a magnetic 



field. For the case of incquivalcnt chains the analysis is 



similar but more involved, and will be discussed else- 
where f34j] . However, we find also in this case indications 
for quasi-long-range order of the pairing order parameter. 

The generalization of the above analysis to AT- leg Hub- 
bard ladders can be performed along the same lines, 
and the presence of a non-trivial plateau when Nn^i = 
n ± m € Z could be predicted. It is important to stress 
however that for the operator responsible for such a 
plateau to be relevant a fair amount of interchain cou- 
pling would be needed. This situation is reminiscent of 
what was found in the study of N— leg spin ladders, where 
the interchain coupling should be strong enough for the 
appearance of non-trivial plateaux [4|. Details of this 
analysis will be presented elsewhere [ 34j . 

The behavior of a 2-leg doped Hubbard ladder in the 
presence of a magnetic field generalizes the 2-leg spin lad- 
der case. Our main result is that the magnetization curve 
should present two plateaux for a non-trivial value of the 
magnetization (see Fig. 1). The plateau at zero magneti- 
zation persists for non zero doping, while the doping de- 
pendent one can be obtained only by keeping the doping 
of the system fixed (as it is the case in most experimen- 
tal settings) . This should be contrasted to the dimerized 
Hubbard chain, in which the zero magnetization plateau 
at half filling is shifted to a non-trivial value when the 
system is doped. More importantly, the doping depen- 
dent plateaux are predicted for a wide range of couplings, 
which should allow for their experimental observation at 
quite low magnetic fields. 
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